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ABSTRACT 

We have investigated the degree to which turbulence in the Very Local Inter- 
stellar Clouds resembles the highly-studied turbulence in the solar corona and the 
solar wind. The turbulence diagnostics for the Local Clouds are the absorption 
line widths measured along 32 lines of sight to nearby stars, y ielding measure- 



ments for 53 absorption components (IRedfield and Linsky II2004I ). We have tested 



whether the Local Cloud turbulence has the following properties of turbulence in 
the solar corona or the solar wind: (a) velocity fluctuations mainly perpendicular 
to the average magnetic field, (b) a temperature anisotropy in the sense that the 
perpendicular temperature is larger than the parallel temperature (or at least en- 
hanced relative to expectation), and (c) an ion temperature which is dependent 
on the ion Larmor radius, in the sense that more massive ions have higher tem- 
peratures. Our analysis of the data does not show compelling evidence for any of 
these properties in Local Cloud turbulence, indicating possible differences with 
heliospheric plasmas. In the case of anisotropy of velocity fluctuations, although 
the expected observational signature is not seen, we cannot exclude the possibility 
of relatively high degrees of anisotropy (anisotropy parameter e ~ 0.50 — 0.70), if 
some other process in the the Local Clouds is causing variations in the turbulent 
line width from one line of sight to another. We briefly consider possible reasons 
for differences between coronal and solar wind turbulence and that in the Local 
Clouds. The apparent absence of anisotropy of the velocity fluctuations and ion 
temperature might be due to randomization of the interstellar magnetic field on 
spatial scales ~ 10 parsecs, but this would not explain the absence of ion mass- 
dependence in the ion temperature. A likely explanation of all the results is the 
greater collisionality, due to ion-neutral collisions, of the partially-ionized Local 
Cloud plasma. 



Subject headings: interstellar medium:magnetic fields 
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Introduction 



Turbulence is believed to be an extremely important phenomenon in many astrophysical 
media, such as the diffuse phases of the interstellar medium, molecular clouds, supernova 
remnants, accretion disks around compact objects, extragalactic radio sources, and the in- 
tracluster media of clusters of galaxies. It is speculated that this turbulence may serve such 
functions as providing additional sources of pressure or heat input, determining transport 
coefficients such as viscosity and resistivity, and accelerating charged particles to high en- 
ergy. Assessment of these possibilities requires knowledge of the properties of the turbulence 
in some detail. To address these possible roles of turbulence, we need to know not only the 
power spectrum, which m ay be nearly i denti cal for turbulent excitations in vastly different 



cases (see, for example, iBayley et al I Il992[ ). but more importantly, the symmetries and 



relationships between fluctuations in plasma velocity, magnetic field, density, etc. 

Such knowledge is almost unattainable for most media of interest to astronomers. In 
many astrophysical plasmas, even the average plasma parameters are incompletely known 
or totally unknown. Information on the properties of turbulence is almost always highly 
indirect in the spatially-averaged, path-integrated measurements available to astronomers. 
Given this, a-priori knowledge of the nature of turbulence from independent sources becomes 
very important. 

Observations of the magnetohydrodynamic (MHD) turbulence which exists in the solar 
wind play a very important role in our understanding of turbulence, and solar wind data sets 
have been the primary data sources in the development and validation of theories of MHD 
turbulence. In this paper, we consider the properties of turbulence in another medium which 
is, perhaps surprisingly, very well diagnosed. This medium is the partially-ionized plasma 
contained in clouds in the Very Local Interstellar Medium (VLISM) within about 15 parsecs 
of the Sun. The properties of these clouds, and the ways in which turbulence in them is 
measured, are discussed in Section 3. The goal of this paper is to investigate the extent to 
which the turbulence in these clouds resembles, or differs from, the turbulence in the solar 
corona and the solar wind. 



A preliminary report on this topic was given in ISpangler et al I (120101) . That paper 



pointe d out that the high spectral resolu tion absorption line measurements of lRedfield and Linsky 
( 1200 ll ) and iRedfield and Linsky I ( 120041 ) could be used to extract properties of turbulence in 
the Local Clouds, and t hat the inferred proper ties could be compared with those of the solar 
corona and solar wind. ISpangler et al (2010) drew p r elimin ary conclusions from examina- 
tion of the results published in IRedfield and Linsky I (120041) . In the present paper, we use 
detailed and quantitative analyses of the data of IRedfield and Linsky I ( 120041 ) to study the 
plasma turbulence in the Local Clouds. 
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2. Properties of Solar Wind Turbulence 

The obvious advantage of solar wind turbulence is that basic plasma physics measure- 
ments of vector magnetic field, plasma flow velocity, density, temperatures, and even electron 
and ion distribution functions can be measured in situ with spacecraft. The fluctuations in 
all of these quantities have been extensively studied, a large literature written, and major 
conclusions reached. A mong the many in fluential articles a nd re views of the subject are 



Bavassano et al ( 



by lGoldstein et al 



982), the m onograph by 



( 119951 ) and iBruno and Carbone I (120051 ) . 



Tu and Marsch I (119951 ). and the review articles 



Another nearby plasma with extensive diagnostics (although not, as yet, in situ mea- 
surements) is the solar corona. Our knowledge of the corona and its turbulence results from 
high spatial resolution images, ultraviolet spectroscopy of numerous transitions, and radio 
propagation measurements. In addition, a sort of "ground truth" of coronal plasma mea- 
surements is provided by spacecraft measurements at heliocentric distances from 0.28 to 1 
astronomical units (AU). The coronal plasma is convected out into space and becomes the 
solar wind. Among the many reviews of the cor o nal p lasm a, two which are par t icular ly 
relevant to the present investigation are ICranmer I (120021 ) and iBird and Edenhofer I (ll990l ). 



A list of the main properties of solar wind and coronal turbulence could be extensive. 
We list four properties which are particularly relevant to the present investigation. 



1. The fluctuations in magnetic field and plasma flow velocity are highly correlated. The 
equations of magnetohydrodynamics couple fluctuations in magnetic field and plasma 
flow velocity. The dimensionless amplitude of the velocity fluctuations, Sv/Va is highly 
correlated with the dimensionless amplitude of the magnetic field fluctuations, Sb/B , 
where 5b and 5v are the root-mean-square (rms) fluctuations in magnetic field and 
flow velocity, re spectively, and Va and Bp a re the Alfven speed and magnitude of the 
magnetic field ( ISpangler and Spitler 112004 ). This property can be used to infer the 
magnitude of fluctuations in one of the quantities (5b or 5v), given a measurement of 
the other. 

2. The fluctuations in magnetic field and velocity are predominantly perpendicular to the 
large scale interplanetary magnetic field. This property can be readily understood if 
the solar wind turbulence is viewed as an ensemble of interacting Alfven waves, or if it is 



descr ibed by the equations of quasi-2D magnetohydrodynamic s (IZank and Matthaeus 



19921) . Observat i ons il lu strating this property a re pr esented in 



Bavassano et al 



Klein et al I (ll993f ) . iBavassano et al I ( 119821 ) and Klein et al I dl993h also show that 



(1982) 



and 

transverse velocity fluctuations are not a universal property of solar wind turbulence. 
They find that the anisotropy of the turbulence (i.e. the excess of perpendicular over 
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parallel fluctuations) decreases with heliocentric distance, and is less pronounced in 
the slow speed solar wind relative to the high speed wind. 

3. A plasma heating process is occurring which preferentially increases the perpendicular 
(to the magnetic field) temperature relative to the parallel temperature. This is a 
hallmark of ion heating by ion cyclotron resonance mechanisms. In the case of the 
corona, the perpendicular temperature T\ exceeds the parallel temperature T\\ by a 
large factor f jCranmer I l2002t iHollweg II2008I ). In the case of the solar wind, typically 
T± < Tj|, but the observed values are much higher than woul d be expected in a n 
expanding solar wind without preferential perpendicular heating (IKasper et alll2009l ). 
This is often expressed as a systematic increase of the fi rst adiabatic invariant with 
heliocentric distance in the solar wind. IChandran I (120101 ) has recently proposed that 
another heating mechanism termed "stochastic acceleration" , due to significant velocity 
fluctuations on scales comparable to the ion cyclotron radius, can also produce the 
observed features of anisotropic temperatures. Chandran's mechanism does not require 
a cyclotron resonance between ions and plane Alfven waves. Regardless of the correct 
theoretical explanation, the aforementioned phenomena are most pronounced for the 
solar corona and the high speed solar wind, i.e. the most collisionless parts of the 
heliosphere. Finally, it is worth noting that in collisionless helio spheric plasmas, the io n 



and electron distribution functions are not Maxwellians (e.g. iTu and Marsch Ill995l ). 
a fact which m ay be of great significance to the thermodynamics of these plasmas 
rtScudder |[l992h . 



4. Different ions, characterized by different Larmor radii, have different temperatures 
Those ions with larger Larmor radii have higher temperatures. This property , which is 



also r eadily understood on the basis of ion cyclotron reso nance mechanisms ( IHollweg 



20081 ) as well as the stochastic acceleration mechanism of IChandran I (120 lOf ). is highly 



pronounced in the solar corona, where the spectro scopically- inferre d temperature of 
OVI is at least 30 times the proton temperature (ICranmer 1 120021 ) . A particularly 
st riking illustrati on of the observed variation of temperature with ion species is given 



m 



Moran I ( 120031 ) . particularly Figure 5 of that paper. 



These properties are prominent in the highly collisionless plasmas of the corona and the high 
speed solar wind, and become les s pronounced in th e denser, slower, and more collisional 
parts of the slow speed solar wind (IKasper et al 1120091 ) . It is of great interest to ask whether 
other astrophysical plasmas also possess turbulence with these properties. A positive answer 
would permit immediate application of the considerable body of knowledge and understand- 
ing of heliospheric turbulence to the numerous media listed at the beginning of this section. 
The specific goal of this paper is a modest step in that direction. We explore the extent to 
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which turbulence in the Local Clouds of the VLISM resembles the turbulence in the solar 
wind and solar corona. 



3. The Clouds of the Very Local Interstellar Medium (VLISM) and the 

Turbulence They Contain 

The Very Local Interstellar Medium (VLISM) is loosely defined as the interstellar 
medium within about 15 parsecs of the Sun. One of the interesting aspects of the VLISM 



is th at it contains about 15 clouds with diameters of a few parsecs (IRedfield and Linsky 



20081 ) . It appears that the Sun is near the interface and region of interaction of two o 



clouds, the Local Interstellar Cloud, or LIC, and the G cl oud (IRedfield and Linsky 



thes e 
2008h . 



Rev iews of the propert ies of these clouds may be found in iFrisch I (120001 ) . IRedfield I (120091 ). 
and iFrisch et al I (120111 ). Most of the information we have about these clouds comes from UV 
and visible wavelength spectroscopy. Absorption lines attributable to these clouds are mea- 
sured along lines of sight to nearby stars with precisely known distances. Properties of these 
clouds are deduced from the Doppler shift, strength, and width of the spectral lines. These 
clouds are plasmas because absorptio n lines of ions as well as n eutral atoms are observed; 
the ionization fraction is about 50 % (IRedfield and Falcon 1 120081 ) . 



Although the information available on these clouds is not as extensive as for the solar 
wind or solar corona, it is sufficient to place the Local Clouds among the best-diagnosed 
astrophysical plasmas. There are several reasons for this state of affairs. First, because the 
absorption lines are measured in the spectra of nearby stars with precisely known distances, 
the spatial extent of the clouds is well determined. Second, the neutral compon ent of the 



cloud s flows into the inner solar system, where it can be measured in situ (e.g. iMoebius 



20091 ). Finally, the heliosphere is embedded in one of these clouds, the LIC cloud, and the 
solar wind interacts with it. The shape and other characteristics of the solar wind interactio n 
provide constraints on the LIC cloud properties (ILallement et al 1 120051 ; lOpher et al 1120091 ). 
The mean plasma properties of the turbulent clouds are given in Table 1 (adapted from 
Redfield and Linsky lliooi IRedfield and Falcon lliooi ). 



Information on turbulence in the Local Clouds is discussed in IRedfield and Linsky 



(120041 ) . Such information is retrievable because the absorption line w idth b can be measured 
for transitions of several atoms or ions. IRedfield and Linsky I (120041 ) fit the line width data 
for each line of sight and Doppler component to the formula 

,2 2 k B T 



ni 



1 



where T is the temperature (assumed the same for all atomic and ion species), k B is Boltz- 
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mann's constant, m is the mass of the atom or ion, and £ is the non-thermal Doppler width 
of the line, attributable to turbulent flows in which all ions and atoms participate. 

To anticipate one of the main points of this paper, a coronal astronomer or solar wind 
physicist would immediately take issue with Equation (1), noting point (4) above that in 
those media, different atoms and ions have different temperatures. In Section 4.4 we will 
investigate the degree to which a single, common temperature characterizes the Local Clouds. 



4. A Study of Local Cloud Turbulence Properties 



In this section, we test whether the four properties of solar wind turbulence stated 
in Section 2 above also characterize the turbulence in the Local Clouds. The data set 



we use are published meas urements of T and £ given in iRedfield and Linsky I (120041 ) and 
Redfield and Linsky I (120011 ). We also have used the l ine width measurements b fo r different 
atoms and ions. These data are shown in Figure 1 of IRedfield and Linsky I (120041 ) . We have 
used the numerical versions of those data, which are available for many lines of sight. 



Data on T and £ are given in IRedfield and Linsky I (120041 ) and IRedfield and Linsky 



( 1200 ll ) for 32 lines of sight, possessin g 53 absorption line compon ents. Data on 50 absorption 
components are given in Table 1 of IRedfield and Linsky I (120041) . Data for the remaining 3 
components are given in Table 5 of IRedfield and Linsky I (1200 ll ). Since some of the 32 lines of 
sight intercept more than one cloud, we have a larger number of absorption components than 
lines of sight. Each of the 53 components provides an independent observational estimate 
of the properties of turbulence in one of the 15 Local Clouds. In the following, we analyse 
these data and investigate the degree to which the Local Cloud turbulence adheres to the 
characteristics listed in Section 2. 



Table 1. Mean Plasma Parameters of Local Interstellar Clouds 



Plasma Parameter Value 

electron density 0.11 cm -3 

neutral density 0.1 cm -3 

temperature 4000-8000 K (typical) 

magnetic field 3-4 /i G (assumed) 
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4.1. Velocity Fluctuations Perpendicular to £> 

We interpret the turbulent line width parameter £ to be a measure of the velocity 
fluctuations in the cloud turbulence. If solar wind turbulence is a good model for the cloud 
turbulence, this turbulence is Alfvenic, and the velocity fluctuations should be perpendicular 
to the large scale interstellar magnetic field B . If this is the case, the measured value of £ 
should vary with position on the sky. 

The reasoning behind this statement is illustrated in Figure 1. In certain directions on 




Fig. 1. — An illustrated argument why the turbulent broadening parameter £ should depend 
on direction on the sky, if the turbulence is Alfvenic with velocity fluctuations _L to the large 
scale interstellar magnetic field B . The heavy arrows indicate the direction of B , and the 
shaded regions indicate the relatively dense Local Clouds, embedded in the rarefied Local 
Cavity. A line of sight along the large-scale field (to star A) will show little turbulent line 
broadening because the fluctuations are perpendicular to the line of sight. A line of sight 
across the field (to star B) will show large turbulent line broadening becaus e the turbulent 
veloci ty fluctuations are aligned with the line of sight. Figure taken from ISpangler et al 



(12010h 



the sky, we are looking across B , and the turbulent fluctuations should be along the line 
of sight (more properly, one of the two fluctuating magnetic field components in a plane 
perpendicular to -Bo will be aligned with the line of sight). For other directions, we are 
looking along B$, and the turbulent velocity fluctuations are mainly transverse to the line 
of sight. In this case, £ should be small. 
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The direction of Bq must be considered an unknown parameter. Although there is infor- 



Minter and Spangler Ill996l ; IVan Eck et al 112011 



Rand and Kulkarni 


1989; 


Rand and Lvne 


1994 



as well as measureme nts of the polarized 



Galactic synchrotron emission ( iBeck et al 1119961 ; lHaverkorn et al Il2004al ). the magnetic field 
models are deduced from measurements on lines of sight which are kiloparsecs in length. 
All analyses of the galactic magnetic field agree that the fluctuating (presumably turbulent) 
component of the galactic magnetic fiel d is comparable to or larger than the systematic com- 
ponent (e.g. iRand and Kulkarni Ill989l ; iMinter and Spangler Ill996l ; lHaverkorn et al 112004a] ; 



Van Eck et al 



201 ll ). By systematic component, we mean a vector field which is describable 



by a relatively simple function of Galactocentric coordinates. 

In a very local sense, the "large scale" magnetic field is almost certainly dominated by 
these turbulent fluctuations. It may reflect the random orientation of the largest eddy in the 
solar neighborhood. In any case, Bq, and the unit vector in the direction of Bo, b, can point 
in any direction in the sky. 

Before proceeding further, it is necess ary to note that there a re two independent, and in- 
compatible estimates of the direction of b. lLallement et al I (120051 ) use the difference between 
the direction of neutral helium flow and that of the largest concentration of neutral hydrogen 

outside the hehopaus jj] to infer that b points in the range I = 205° — 240°, b = —60° 38°. 

Gurnett et al I (120061 ) report measurements of the direction to sources of low frequency radio 
emission, which are assumed to be generated in the heliosheath @. They assume that this 
radio emission is generated at points on the heliopause which are perpendicular to B . They 
do not retrieve the vector b, but report the angle between b and the direction to t he ecliptic 
pole. They interpret their results as being consistent with lLallement et al I (120051 ). 



A mo del- dependent estimate of b has also been presented by lOpher et al I (120091 ). lOpher et al 



(120091 ) use Voyager 1 measurements of the plasma flow direction in the heliosheath, together 



I = 1 0° - 20°, b = 28° - 38°. 



( 


2005 


) and 


Opher et al 



(120091 ) appear to be significantly different. It is worth emphasizing that all of the afore- 
mentioned techniques are model-dependent in that they adopt physical assumptions about 
processes in the outer heliosphere, or use MHD simulations of the heliosphere to relate the 
actual measured quantity to the properties of the solar wind-ISM interaction, including the 



lr The heliopause is the contact discontinuity between the shocked solar wind and the interstellar medium. 

2 The heliosheath is the region of shocked solar wind between the solar wind termination shock and the 
heliopause. 
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direction of b. 



To carry out an analysis suggested by Figure 1, we need a model for the form of the 
turbulent line width £ in the case of anisotropic, Alfvenic turbulence in the Local Clouds. 
The derivation of such an expression is given in the Appendix. We assume the turbulence is 
characterized by a root-mean-squared amplitude of the perpendicular velocity fluctuations 
V±, and an anisotropy factor e. When e is zero, the turbulence is isotropic, and if e = 1 the 
turbulent motions lie completely in the plane perpendicular to B , with no motion in the 
direction of the field. In this case, (see Appendix for derivation) 



e_ 



V? 



1 — e(sin b sin f3 + cos Al cos b cos j3)* 



(2) 



In addition to variables already defined, the galactic coordinates of the direction of the line 
of sight are (/, b) and the galactic coordinates defining the direction of the local magnetic 
field are (A,/3). The angle Al = A — I. The component of the gas velocity along the line of 
sight is given by vl- Equation (2) is useful for analyses of the sort to be described shortly 
but a more intuitive expression is 



e_ 



< v L > 
V? 



1 - e cos 2 A 



(3) 



where A is the angle between the line of sight and the local field. 



We cannot, in a straightforward way, test whether actual data adhere to the relationship 
in Equation (3) without knowledge of the direction of the local magnetic field, indicated 
by the unit vector b or the angles (A, 0). We therefore adopted the following, brute-force 
approach, which was under taken without re f erence to t he a-priori es t imate s of the magnetic 
field direction proposed by lLallement et al I (120051 ) and lOpher et al I (120091 ) . 



1. We chose 361 candidate directions for b, each characterized by values of (A, 0). These 
candidate directions were spaced 10 degrees apart in galactic latitude and longitude, 
and completely covered one hemisphere of the sky. Given the nature of the anisotropy 
sought, one hemisphere is adequate for complete coverage. 

2. For each candidate direction, the angle A could be calculated (u sing Equation (2)) for 



each l ine of sight for which estima t es of T and £ are available from lRedfield and Linsky 



( 20041 ) and iRedfield and Linsky I ( 20011 ) . We then made a plot of ^ versus cos A for 
this candidate direction of b. 



3. 



Each plot was examined to see if a relationship of the form given in Equation (3) could 
be discerned. 
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No compelling cases for such a relationship were found. A set of cases in which something 
like the expected relationship seemed to be present in the data (often referred to by self- 
deluded individuals as "tantalizing") were collected for further scrutiny. In these cases, 
there seemed to be a larger average value of £ 2 for smaller values of | cosA| than for larger 
values. The total number of such cases was 10, and these cases were roughly clustered in 
the direction A ~ 40°, f3 ~ 60°. The case for (A = 40°, = 60°) is shown in Figure 2. We 

O 



CM 



* O 




0.2 0.4 0.6 0.8 1 

COS(A) 

Fig. 2. — Plot of £ 2 versus | cos A\, when the local interstellar magnetic field is assumed to 
point in the direction A = 40°, (5 = 60°. The plot shows measurements of £ 2 extending to 
larger values for cos A < 0.40, as expected for transverse velocity fluctuations. The statistical 
significance of this difference is discussed in Section 4.1.1. The dashed line shows the model 
given by Equation (3) with V]_ = 8.00 km 2 /sec 2 , and e = 0.70. 

do not claim this result as a detection of anisotropy in the Local Cloud turbulence, but 
it does give an indication of the data quality in one of the best cases for anisotropy. An 
obvious feature in Figure 2 is the presence of two data points at | cos A\ ~ 0.8 — 0.9 with 
possibly anomalous values of £ 2 . These points correspond to the +2.6 km/sec radial velocity 
component for Alkaid (HD 120315) and the +13.9 km/sec radial velocity component for ( 
Doradi (HD 33262). Obviously, any hint of a systematic dependence of £ 2 on A disappears 
if these are valid points. However, it appears that they lie well outside the distribution of £ 2 
values for stars with similar values of cos A. There is basis for suspecting that the line widths 
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might be affected by biends of two or more components, thus inflating a single component 
fit to the line. For this reason, Alkaid and £ Dor may be provisionally considered outliers. 
In the analyses which follow, we have considered the complete data set of 53 absorption 
measurements, as well as an edited subset in which Alkaid and £ Doradi are removed. 



4-1.1. A Search for Weaker Anisotropy and Quantitative Limits to the Anisotropy 

Figure 2 does not present a strong case for anisotropy, defined as close adherence of the 
data to the expression given in Equation (3). However, it is possible that an anisotropy of the 
sort we are seeking is present, but obscured by star-to-star variations of another, unknown 
nature. To detect anisotropy in this case, it is necessary to average measurements for several 
stars. Furthermore, we need a means of extracting from the data a quantitative upper limit 
to the anisotropy parameter e. 

A simple way of addressing both of these points is to average the data over intervals 
in cos A For these purposes, we consider the star-to-star variations as noise superposed 
to a true signal of the form in Equation (3). If anisotropy is present, the average value of 
fpz for all lines of sight with < cos A < 0.3 will be larger than for all lines of sight with 
0.7 < cos A < 1.0. Furthermore, the ratio of the mean values for £ 2 in the two ranges of cos A 
is a measure of, or upper limit to, the anisotropy. The following analysis was undertaken. 



For each of the 10 lines of sight for which there was some suggestion of anisotropy, as 
in the case of Figure 2, we computed a list of £ 2 versus cos A Once again, we point 
out that an assumed, candidate direction for the local interstellar magnetic field is 
necessary to calculate the angle A 

The mean value of £ 2 was calculated for all measurements in two ranges in cos A 
< cos A < 0.3, and 0.7 < cos A < 1.0, as well as an estimate of the error in the mean. 
The error in the mean of £ 2 in the intervals was calculated as follows. We used the 
measure d mean and standard deviatio n of £ values (2.24 and 1.03 km/sec respectively) 



given by iRedfield and Linsky I (120041 ). in Figure 2 of that paper. These values were 
used to compute the mean and standard deviation of the quantity £ 2 . Finally, the 
standard deviation of the mean of the quantity £ 2 for a sample of 16 measurements (30 
% of 53 data points) was calculated. 

A ratio R was calculated in which the numerator was £ 2 =< £ 2 >i in the first interval, 
and the denominator was £ 2 =< £ 2 > 2 in the second interval. The error in this ratio 
was calculated in the standard way, using the standard deviation of the mean of £ 2 



- 12 - 



described in item # 2 above. Our value for the standard deviation in i?(0.3), used for 
all directions, was 0.33. 

4. Steps # 2 and # 3 were repeated for the broader intervals of < cos A < 0.5 and 
0.5 < cos A < 1.0 (a = 0.5). The associated error in R(a = 0.5) was taken to be 0.26. 

The reason for carrying out the ratio analysis for two values of the interval width a, 
a = 0.3 and a = 0.5, is as follows. An analysis of this sort has competing demands on the 
value adopted for the width of the averaging interval, a. The smaller the value of a, the 
greater will be the contrast between the mean values of < £ 2 > for the two intervals. On the 
other hand, a larger value of a results in more stars and absorption components in each bin, 
and thus a statistically more stable value of < £ 2 >. We carried out an analysis for both 
a = 0.3, to try and get the largest possible contrast with a significant number of data points 
contributing to the average, as well as a = 0.5, which has lower contrast but includes all the 
data in the sample. 

The mean value of £ 2 in each of the intervals, and their ratio, is easy to calculate from 
the data. To relate this ratio to the anisotropy factor e requires use of Equation (3) for the 
expected relationship £ 2 (A). Given Equation (3), the mean value < £ 2 >iin the first interval 
< cos A < a is given by 

T/2 ra 

< e > 1= — / (1 - ex 2 )dx (4) 
a Jo 

and similarly with the second interval, so the ratio of the two means, R(a) is given by 
m \ _ < e >i Jp a (l - ex 2 )dx 1 - ea 2 /3 

Ria) = = — r = — — — = Ria.e) 5 

K J <e>2 fl a (l - tx^)dx 1 - e(a 2 - 3a + 3)/3 1 ' { J 

A plot of R(a, e) for a = 0.3 and a = 0.5 is shown in Figure 3. 

Values of R(a = 0.3) and R(a = 0.5) and associated anisotropy indices e were calculated 
for the 10 candidate directions for b described in Section 4.1. Data from all 53 lines of sight 
were used. In only one case, for A = 10°, /3 = 40°, did R(a = 0.3) exceed unity by more than 
twice the adopted error (R = 1.71 ± 0.33), and even in this case the value of R(a = 0.5) was 
consistent with unity. For the other nine directions, both i?(0.3) and _R(0.5) were consistent 
with unity (i.e. within 2 standard deviations of unity). Of the twenty calculated quantities 
(R(0.3) and R(0.5) for 10 directions), half had R < 1, which is inconsistent with a velocity 
anisotropy of the sort we are seeking, but is consistent with random variations in the case 
when £ 2 is isotropic. 

The results of this analysis may be summarized as follows. If the data for Alkaid and ( 
Doradi are valid values, then there is no candidate magnetic field direction with a statistically 
significant value of R > 1, and corresponding value of e different from 0. 
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0.2 0.4 0.6 0.8 1.0 



Fig. 3. — Plot of the ratio R(a,e) for a = 0.3 (solid line) and a = 0.5 (dashed line). Given 
an interval for averaging data (a = 0.3, 0.5), the measured ratio of the mean values of £ 2 in 
the intervals < cos A < a and 1 — a < cos A < 1 corresponds to the ordinate, and the 
inferred value of the anisotropy factor e is the abscissa. 



We now repeat this analysis, but assuming that Alkaid and ( Doradi are outliers which 
may be excluded from the sample. The results are presented in Table 2. The first four 
columns of this table contain, respectively, the galactic latitude and longitude of the can- 
didate field direction, the measured value of R with a = 0.3 and associated error, and the 
value of the anisotropy parameter e consistent with R and its error, obtained from Equation 
(5). The final two columns give the value of R with a = 0.5 and the value of the anisotropy 
corresponding to this value of R. The first ten rows correspond to the directions selected 
from our visual examination of plots similar to Figu re 2. The candidate field directions in the 



bottom two rows correspond to those proposed by lOpher et al I (120091 ) and lLallement et al 



(120051 ). respectively, and are discussed further in Section 4.3. For the moment, we restrict 
consideration to the first 10 rows of Table 2, which were directions (A, (3) chosen by us for 
closer examination. 

Table 2 shows a number of candidate directions for b with marginally significant evidence 
for anisotropy. Directions such as (A = 40°, (3 = 40°), (A = 40°, (3 = 60°) (the case shown 
in Figure 2), (A = 60°, /3 = 50°), and (A = 90°, (3 = 50°) have values of R{a = 0.3) which 
exceed unity by 2a or more. A plot of £ 2 vs | cos^4| for (A = 90°, (3 = 50°) is shown in Figure 
4, in the same format as Figure 2. 

In some of these cases both R(a = 0.3) and R(a = 0.5) exceed unity by about 2a or 
greater. Furthermore, the inferred values for e for the 2 binning intervals are in agreement, 
within the errors. Although the results of this analysis do not present strong evidence for 
anisotropy of the Local Cloud turbulence (our best cases are, after all, shown in Figures 2 and 
4), they are not inconsistent with b pointing in the direction (A = 40° ± 20°, (3 = 50° ± 20°), 
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Fig. 4. — Plot of £ 2 versus | cosA| for A = 90°, /3 = 50°. The plot is in the same format as 
Figure 2. The dashed line shows the model given by Equation (3) with V]_ = 6.0 km 2 /sec 2 , 
and e = 0.75. 

and an anisotropy parameter e = 0.5 — 0.7. Before leaving this section, two points should be 
emphasized. First, the modest indications of anisotropy in Table 2 are completely dependent 
on excluding the measurements of Alkaid and £ Dor, which so prominently depart from the 
model curves in Figures 2 and 4. Second, the measurements of £ 2 vs. | cos A\ do not adhere 
closely to the relationship given by Equation (3), but show a dispersion about that curve 
which is larger than the measurement error. If anisotropy is present in these data, there 
must be another, unnamed physical process responsible for variation in £ from one line of 
sight to another. 



4-1.2. The Value of e for the Solar Wind 



Having presented our results on upper limits to the anisotropy factor for turbulence 
in the Local Clouds, we now consider the corre sponding quantity in the sol ar wind, which 
has the benefit of direct, in-situ measurements. ISpangler and Spitler I ( 120041 ) resolved solar 
wind magnetic field fluctuations into components parallel and perpendicular to the large- 
scale interplanetary magnetic field. The data came from the magnetometer of the WIND 
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spacecraft § at a heliocentric distance of about 1 AU. Spangler and Spitler ( 2004 ) used 66 
intervals of one hour duration during slow solar wind conditions, and report their results in 
terms of modulation indices m B \\ and m^i of fluctuations parallel and perpendicular to the 
mean field 0, 



m B \\ 
m B± 



Sb\\ 
Sbi 



B, 



o 



(6) 
(7) 



where 5b» and 5b± are the rms values of the fluctuations in th e magnetic field components 



parallel and perpendicular, respectively, to the mean field B Q . iSpangler and Spitler I (120041 ) 
report mean values for m_Bii and m B ± of 0.0321 and 0.112, respectively. The means are for the 
distrib ution of values measured in the 66 data intervals. As noted in ISpangler and Spitler 
(120041 ). rriB± should be larger than m B \\ because it possesses contributions from two turbu- 
lent field components rather t han just one. The degree of anisotropy can be determined by 
comparing rriB\\ to m B ±/\/2. ISpangler and Spitler I (120041 ) found that the turbulent fluctu- 
ations in their study were Alfvenic, in the sense that |^ = Jp, so the measured magnetic 
field modulation indices may be considered proxies for modulation indices of the velocity 
fluctuations. With this assumption, we have 



V2m 



B\\ 



1- ^ 



m B ±_ 

2 



V \\ 



2m 2 m 



711 



(9) 



B± 



where in Equations (8) and (9) we make the connection between the root-mean-square ve- 
locity fluctuations and the velocity scales V\\ and V± of the fluctuation distribution funct ion 
in the Appendix. Using the values of m B \\ and m B ^ from ISpangler and Spitler I (120041 ) in 
Equation (9), we have e = 0.84. 



The anisotropy of solar wind fluct uations had been co n sidered prior to ISpangler and Spitler 

d2004h bv lBavassano et all Jl982f ) and lKlein et~aTl dl993h . iBavassano et al I dl982h also stud- 
ied magnetic field fluctuat ions, in condition s of hi gh speed solar wind at several heliocentric 

( 119821 ) (see data in their Figure 2) yield values 



distances. The results of Bavassano et al 



3 WIND is one of the spacecraft which comprises the International Solar- Terrestrial Physics (ISTP) pro- 
gram. 



Spangler and Spitler I ([2004) used the variables e^ii and cb± for the modulation indices, but we do not 



retain this notation so as to avoid confusion with our anisotropy parameter e. 
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for e in the range 0.8 - 0.9 and greater, i.e. very similar to that quoted above. It should 
also be noted that a highly anisotropic and Alfvenic nature is a characteristic of turbu- 
lence in the inner solar syst em, that mi g ht no t be valid throughout interplanetary space. 



Bavassano et al I (119821 ) and Klein et al I (119931 ) show that properties of solar wind turbu- 



lence depend on heliocentric distanc e, specifically that the degree of anisotropy decreases 
with increasing heliocentric distance. iKlein et al I (119931 ) further report that the anisotropy 
is less in slow speed than high speed solar wind. The upper limit we can place to anisotropy 
of velocity fluctuations in the Local Clouds is less than, though comparable to, that of solar 
wind turbulence at a heliocentric distance of 1 AU. 



4.2. Anisotropy in the Ion Temperature 



As discussed in Section 3, the analysis of line widths by iRedfield and Linsky I (120041 ) 
also yields the ion temperature T. Strictly speaking, this is a line-of-sight temperature; it is 
a measure of the line-of-sight component of thermal motion of atoms and ions. If the Local 
Cloud turbulence is similar to heliospheric turbulence, the ion temperature T might also 
depend on the angle A between the line-of-sight and the local interstellar magnetic field, as 
a consequence of T± ^ Tji. In regions of the solar wind in which ion cyclotron resonance 



heati n g appears to be acti ve, and in particular in the solar corona, T± ^> Tji ( ICranmer 



20021 1 . iKasper et al I (120091 ) have made an extensive study of the perpendicular-to-parallel 
temperature ratio R p = ^ for protons in the solar wind. Although it is not the case that 
R p > 1 in all cases, it is tr ue that Tj is larger than it would be in an adiabatically-expanding 
solar wind. Furthermore, IKasper et al I (120091 ) identify boundaries in a (R p , 0) plane, where 
(3 is the conventional plasma physics quantity of the ratio of thermal to magnetic pressure. 
These boundaries are clearly seen in the distribution of measured values of R p and j3, and 
correspond to instability thresholds for generation of plasma waves. Apparently these plasma 
waves heat the protons in a way which keeps them "in bounds" in a subset of the (R p , (3) 
plane. The important point here is that in the corona and solar wind, the perpendicular 
and parallel temperatures are usually not equal, and Tj_ > Tji where ion cyclotron resonant 
effects are important. 

To test for temperature anisotropy in the Local Clouds, an analysis similar to that 
of Section 4.1 was undertaken. Plots of T as a function of cos A were made for all 361 
candidate directions for b. These plots were visually examined for indicators that the data 
were organized according to a relation like that in Equation (3). In the case of temperature 
anisotropy, the line-of-sight temperature should obey a relationship like Equation (3). No 
case of a convincing temperature anisotropy was found. Figure 5 shows the results for the 



-re- 



direction (A = 10°, = 40°). 



o 




COS(A) 

Fig. 5. — The temperature as a function of | cos^l in the case of (A = 10°, = 40°). There 
is no obvious sign of temperature anisotropy, in the sense that T± > T\\. The dashed line 
represents a relationship similar to Equation (3), in which T± = 8A25K, and e = 0.45. 



An analysis similar to that of Section 4.1.1. was undertaken, in which a parameter 
Rt was calculated for the mean temperature in two distinct intervals of width a. Values 
for the ratio Rt and associated errors were calculated for the same candidate directions as 
in Table 2. The stars Alkaid and £ Doradi were not excluded from this analysis because 
they are not anomalous as regards temperature. In none of these directions was there a 
convincing case for temperature anisotropy. In only one of the cases examined, that of 
(A = 10°, = 40°) shown in Figure 5, did both Rt(cl = 0.3) and Rx(a = 0.5) exceed unity 
by an amount that was more than twice the adopted error. For this direction, we calculate 
Rt(cl = 0.3) = 1.45 ± 0.08 and Rt((i = 0.5) = 1.2 3 ± 0.06. The errors were ca lculated from 
the dispersion in the measurements of T given in iRedfield and Linsky I (120041 ). in a manner 
similar to that used in Section 4.1.1. These values for Rt would correspond to anisotropy 
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parameters e = 0.34 — 0.44. We do not claim this direction clS db CclSG for temperature 
anisotropy because the data shown in Figure 5 do not show clear adherence to Equation (3). 

For the other directions, the data indicate factors closer to unity, and anisotropy 
factors closer to zero. We thus find no evidence in the data for an anisotropy in the sense 
T± > Xj|, although (as illustrated by the discussion in the previous paragraph) we cannot 
exclude the possibility that e < 0.40 could be present, but hidden by random variations in 
T from one line of sight to another. An upper limit to the temperature anisotropy e < 0.40 
corresponds to ^ < 1.67, which is considerably less than that reported for the solar corona, 



and also less than many cases reported in the solar wind (see Figure 1 of iKasper et al 



20091 ) . 



4.3. Analysis of Data for Select Candidate Magnetic Field Directions 



The analysis of Section 4.1 was done in an unbiased fashion, i.e. with no a-priori 
estimate of the local direction of the interstellar magnetic field. No direction examined had 
a compelling case for anisotropy of the turbulent amplitude £ or the ion temperature T. With 
this analysis co mpleted, we then re-ex amin ed the data for "pre ferred" candidate directions 
b advocated by lLallement et al I (120051 ) and lOpher et al I ( 120091 ). as discussed in Section 4.1 
above. 



Lallement et al I ( 120051 ) propose a direction of the local interstellar magnetic field of 
205° < A < 240°, —60° < j3 < —38°. Taking the means for each coordinates, we have a 
candidate field direction for lLallement et al I H of (A = 222.5°, p = -49°). Since -b 
serves equally well as a direction for anisotropy, we have (A = 42.5°, (3 = 49°) as a candidate 
direction for the local field. It is interesting that this direction is very close to the set 
of directions chosen as the best candidates from the unbiased analysis of Section 4.1, and 
illustrated in Figure 2. Although the weak, if not nonexistent, evidence of anisotropy in 
Figure 2 precludes any further claims, the coincidence of the set of directions cho s en for 
closer examination and the proposed direction of the local field of lLallement et al I (120051 ) 
could motivate future investigations with more lines of sight. 



The direction for b proposed by lOpher et al I (120091 1 . (A ~ 15°, ~ 33°) did not emerge 
from our unbiased analysis of the £ data as one of the directions for closer examination. 
The values of R(a = 0.3 ) and R(a = 0.5) for a direction in the range of possible directions 
chosen by lOpher et al I (120091 ). (A = 17.7°, = 34.1°) are given in Table 2. The value 
of R(a = 0.3) = 1.46 ± 0.33 is not statistically significant, but is not inconsistent with an 
anisotropy of the velocity fluctuations < 0.4—0.5. Interestingly, the anisotropy for R(a = 0.5) 
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is larger, but not in agreement with the results for the smaller binning interval. 



4.4. A Test for Ion Cyclotron Resonance Heating in the Local Clouds 

As noted in Section 2, in the solar corona and solar wind there is not a single temperature 
which is valid for all ions, as written in Equation (1). Indeed, the temperature increases for 
ions with larger Larmor radii. The reason for adopting Equation (1) in application to the 
Local Clouds is the simple fact that it yields entirely satisfactory fits to t he spectral line width 



data f or lines from as many as 8 different atoms and ions (see Figure 1 of lRedfield and Linsky 
J2004K 



The analysis of this section will be in the nature of establishing an upper limit to the ion 
mass dependence of the ion temperature in the Local Clouds. If cyclotron resonant heating 
is occurring, one would expect a modification of Equation (1). A plausible candidate form 
is 

2 2ksTo ( m \ 



m 



\m J 



(10) 



where T and m are the temperature and mass of the lightest atom or ion analysed, and m 
is the mass of the more massive atom or ion. This equation essentially says that the atomic 
or ionic temperature T{m) = To(m/mo) d . The form chosen for Equation (10) is relevant 
because, in the solar coronal case, the heating has been shown to be more pronounced than 
"mass proportional heating" (jCranmer 1120021 ). which corresponds to d > 1. 



A fit of Equation (10) to the data introduces three model paramet e rs (Tn ,£,cO, rather 
than the two parameters of Equation (1) utilized by lRedfield and Linsky I (12004 ) . This means 
that there is a broader basin of acceptability in a x 2 sense. 

The following analysis was undertaken. 



We selected data from all lines of sight and cloud components which possessed 7 or 8 
transitions, including the deuterium (important for determining To) and iron (impor- 
tant for determining £) line measurements. These data consisted of the measured line 
widths b and associated errors. We had 11 such absorption components for analysis. 

A least-squares fit of Equation (10) was made to the data, and the range of parameters 
T ,£, and d determined which allowed an acceptable value for the reduced xl- We 
chose a value of X 2 = 2.21, corresp onding to a 5 % probability of occurrence for 5 
degrees of freedom (IBevington Ill969l ). as the limit of an acceptable fit. 
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3. Upper limits to d were chosen which corresponded to the maximum value of d consistent 
with the xl — 2.21 acceptability criterion, with no imposed constraints on T Q and £. 



4-4- 1- Results 



The results of this analysis are given in Table 3. Column 1 gives the star name and 
column 2 gives the absorption component, identified by its velocity. Columns 3 and 4 give 
the temperature and turbulent velocity parameter £ from the 2 pa rameter fit, Equation 



(1). T he numbers in these columns are taken directly from Table 1 of iRedfield and Linsky 



(120041 ). and are reproduced here for comparison with the parameters of the model given by 
Equation (10). The errors in columns 3 and 4 are approximations to those published by 
Redfield and Linsky I (120041 ); those authors allowed for different errors above and below the 
mean value. Columns 5,6, and 7 give the values of T , £, and d for the limiting acceptable 
3-parameter model as defined in point # 3 of the previous section. 

Of the 11 absorption components with 7 or 8 measured transitions, 9 gave acceptable 
fits with the temperature model of Equation (10) and d > 0. For two of the 11 components, 
the minimally-acceptable fit of Equation (10) yielded a physically implausible model for b(m) 
in which b increased with increasing m. We will discuss those cases below. 

For the remaining 9 line of sight/absorption components, a minimally-acceptable fit to 
the data was possible with values of d ranging from 0.12 (Capella) to 0.74 (-20.5 km/sec 
component of i Cap). It is worth emphasizing again that these values represent the max- 
imum value of d which is statistically acceptable as defined above, and with little or no 
constraints placed on the values of T and £. Any larger value of d is incompatible with the 
data. An illustration of a fit of Equation (10) to one of our data sets is shown in Figure 
6. This figure also shows the b est-fitting 2 parameter model, with the parameters reported 
by IRedfield and Linsky I (120041 ) . We do not show the model with the maximally acceptable 
value of d — 0.69, but instead Equation (10) with a slightly smaller value of d — 0.60 (and 
associated parameters T and £) that provides a better fit to the data. In all 9 cases, the 
fit of Equation ( 10) introduced lower values of £ than the two-parameter fits of Equation 
(1) published in IRedfield and Linsky I (120041). F or ex ample in the case of the 8.6 km/sec 
component of G191-B2B, IRedfield and Linsky I (I2004J ) report £ = 3.27 km/sec, while the 
limiting acceptable fit of Equation (10) has £ = km/sec and d = 0.69. In fact, as may be 
seen in Table 3, most of the limiting acceptable fits of Equation (10) to the data have £ = 0. 

The reason for this anticorrelation of £ and d is clear. The quantity £ is determined 
by the degree to which b(m) asymptotes to a constant level as m becomes large. However, 
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Table 2. Averaging and Anisotropy Analysis- Alkaid and £ Doradi Excluded 
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P 
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1.96 ±0.33 


0. 
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.09 


1.24 ±0.26 
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13 
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Table 3. Limits on Mass Dependence of Ion Temperature 



Star 


Comp (km/sec) 


T(K) 


£ (km/sec) 


2b(K) 


£ (km/sec) 
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1.68 ±0.39 


6450 


1.30 


0.12 
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0.0 
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0.0 


0.61 
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5.47 ±0.41 


4296 


0.41 


1.00 
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3.93 ±0.22 


1000 
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Atomic Mass (amu) 

Fig. 6. — Spectral line width data b(m) for the 8.6 km/sec absorption component in the 
spectrum of G191-B2B. Plotted is the measured spectral line width b versus the atomic or 
ionic mass in atomic mass units. Solid symbols represen t ions, open symbols are ne utral 



atoms. The solid curve is the best fit of Equation (1) from iRedfield and Linsky I (120041 ); the 
parameters of the model are T = 4400-fT and £ = 3.27 km/sec. The dashed curve represents 
Equation (10) with T = 4133K, £ = 1.63 km/sec, and d = 0.60. The reduced chisquared 
xt for the solid curve (d = 0) is 0.64. The dashed curve, corresponding to d — 0.60, has 
xt = 1.82, which is larger, but still statistically acceptable by our criteria. 



this behavior can also be produced (over a limited range in m) by a lower value of £, if 
compensated by a higher thermal width for more massive ions (the meaning of the (m/mo) 
term) . 

The conclusion from this analysis is that the data do not exclude a modest dependence of 
the ion temperature on mass, with an associated drop in the turbulence level £. This result 
is somewhat ironic in that ion-mass- dependent temperature was sought as a consequence 
of turbulence or wave-particle interaction. However, if mass dependence of temperature is 
present, the data require a reduced level of turbulence in the Local Clouds. 

Two points should be emphasized in this context. First, there is no statistical argument 
for choosing Equation (10) over the simpler Equation (1) as a temperature model; both 
provide statistically-acceptable representations of the data, and in almost all cases the best 
fit has d — 0. Second, with the exception of the two lines of sight with poor fits, the degree 
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of mass-dependence of the temperature which is compatible with the data is less than the 
case of mass-proportional temperature (d = 1) that is the benchmark for the corona. 

To conclude this section, we briefly discuss the two cases in which a plausible model 
for T{m) according to Equation (10) was not obtained. These cases are the 14 km/sec 
absorption component for C Dor, and the 1.7 km/sec component of v Peg. In the case of 
£ Dor, the value of £ from iRedfield and Linsky I (120041 ) was particular ly high (5.47 km/sec. 



noted above) and in the case off Peg, the temperature was low (1700K. IRedfield and Linsky 



(120041 ) ). In both cases, therefore the measured variation of b with m was relatively small, so 



fit trends with -r- > or -i— < were equally compatible with the data. 

dm dm ~l j t- 



4.5. Additional Remarks on Heating by Plasma Waves 



The previous section dealt with a search for mechanisms that produce more effective 
heating of ions with lower cyclotron frequencies. Two mechanisms have been proposed 
which can explain this, both of which involve the interaction of ions with plasma waves and 
turb ulence. These a re the ion-cyclotron resonance interaction discussed bv lCranmer I (12 002') 
and iHollweg I ( 120081 ) inter aliis, and the stochastic acceleration mechanism of IChandran 
( 120101 ). However, our results on spectral line widths as a function of ionic or atomic mass 
place more general constraints on plasma physics processes operative in the Local Clouds. 
Such processes, like magnetic reconnection, interaction with plasma waves, response to large- 
scale electric fields, etc, act directly only on ions since ions respond to e l ectric and magnetic 
fields, while neutral atoms do not. The data of IRedfield and Linsky I ( 120041 ) include both 
neutral atoms and ions. The data in Figure 6 include measurements for lines of CII, Mgll, 
AM, Sill, and Fell, which are ions (solid symbols), as well as lines of the neutral atoms DI, 
NI, and OI (open symbols). In Figure 6, the two-parameter curve fit s all transitions with 
the sa me temperature. The more extensive set of results in Figure 1 of IRedfield and Linsky 
( 120041 ) yields the same conclusion; a single temperature provides a satisfactory description 
of the data for both ions and neutral atoms. This result is obtained for many lines of sight 
throughout the sky. This indicates that whatever process heats the Local Cloud plasmas to 
temperatures of order 7000K does so via a process which heats ions and neutrals equally. 
Alternatively, if heating power is input preferentially or exclusively to ions via some plasma 
physics process, that power is efficiently shared with neutrals. 



5. Discussion 



In this section, we consider the observational findings of Section 4, and speculate on 
their consequences for our understanding of the turbulence in the Local Clouds. 



5.1. Anisotropy of Velocity Fluctuations 

Our analysis shows no compelling evidence that the turbulent velocity width £ depends 
on direction in the sky, in a way which would be consistent with velocity fluctuations pre- 
dominantly perpendicular to a large scale magnetic field (see cartoon in Figure 1). The fact 
that the few directions chosen for clos er inspection have b in approximately the same direc- 
tion as that independently inferred by lLallement et al I J2005h is an interesting curiosity and 



probably motivates future examinations with more data, but does not change the negative 
nature of the current conclusion. We feel that the isotropy of £ is significant, and con- 
tains information on the nature of turbulence in the Local Clouds. There are three possible 
explanations for this result, which are not mutally exclusive. 

• The simplest explanation is that turbulence in the Local Clouds is not Alfvenic, so that 
the velocity fluctuations are not primarily perpendicular to the large scale field in the 
solar neighborhood. The velocity fluctuations parallel to the large scale field are then 
comparable to those in the two perpendicular directions. The mechanisms responsible 
for decay of turbulent anisotropy in the remote solar wind (cited in Sections 2 and 
4.1.2) should be studied as being potentially relevant to the Local Clouds. 

• A related sugge stion is that the fluctuations are transverse on scales much smaller than 



the outer scale (ISpangler 1119991 ). but that the measurements of £ are dominated by the 
largest amplitude fluctuations on the outer scale. These outer scale fluctuations are of 
sufficient amplitude to nearly totally randomize the average field, and thus eliminate 
the simple observational signature sought for in Figure 1. 

The third explanation is closely related to the second, and would say that turbulent 
fluctuations nearly completely randomize the local Galactic magnetic on a scale com- 
parable to the separation between the Local Clouds. According to this viewpoint, 
Figure 1 would refer only to the average magnetic field, and not to the true value at 
a specific cloud (or even within a cloud). In terms of the mathematical vocabulary 
introduced in Section 4.1, the unit vector b would be randomized on scales comparable 
to the spacing between clouds, if not within individual clouds. 
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The last two explanations would seem to be the most plausible and consistent with our 
understanding of turbulence. It is, after all, the dominant interpretation of why Galactic 
Faraday rotation meas urements give such an ambiguous indication of a l arge scale Galac- 



tic magnetic field (e .g. iRand and Kulkarni I Il989l ; iRand and Lyne I Il994j ; iBeck et al 1 11996 



Van Eck et al Il201ll ). However, in the present context, there are arguments against their 



applicability. 

Acceptance of point # 3 would require the outer scale of the inte rstellar turbulence (at 



l east i n the vicinity of the Sun) to be of the order of a few parsecs. iMinter and Spangler 



( 119961 ) estimate the outer scale of the three dimensional turbulence to be about 4 parsecs, 
but suggest that there is a two dimensional component with an outer scale of order 100 
parsecs. Since the rms ampli tude of the magnetic field flu ctuations is dominated by the two 
dimensional component, the IMinter and Spangler I (119961 ) estimate of the turbulent Galac- 
tic field suggests that the correlation length in b is substantially larger than the distances 
between the stars used in this study. A rigorous analysis has not been done however. 

Estimates of the outer scale in interstellar turbulence whic h may lead to a differe nt 



concl usion have been made by Haverkorn and her co-workers (e.g. lHaverkorn et al 



20061 ); a summary of Haverkorn's results on the outer scale is given in ISpangler et all (boiof ). 



2004a 



b, 



Haverkorn has presented results using both Faraday rotation of background radio sources 
as well as the statistics of polari zation fluctuations in the Galactic polarized synchrotron 
radiation. lHaverkorn et al I (120061 ) claim evidence for a difference in the properties of fluctu- 
ations in the spiral arms of the Galaxy relative to the interarm regions. For the case of the 
fluctuations in the spiral arms, their estimates of the outer scale range from 2-17 parsecs. 
There was no effort in this work to distingu i sh bet ween 2D and 3D turbulence, with different 
outer scales. As noted in ISpangler et al I (120101 ) . an outer scale as small as 2-17 parsecs 
(particularly the low end of that range) would probably randomize the local magnetic field 
sufficiently to eliminate anisotropy of £ on the sky. These considerations provide another 
reason for improved determination of the turbulent outer s cale in the interstellar me dium . 
It s hould be kept in mind that t he observational results of ( Haverkorn et al 2004al lbl. 12006 ) 
and IMinter and Spangler I (119961 ) most probably refer to the Warm Ionized Medium (WIM) 
component of the interstellar medium, and that their results are not obviously transferable 
to the turbulence in the Local Clouds. 

A second argument against points #2 and #3 is that we have information on the 
dimensionless amplitude of the Local Cloud turbulence, an d the amplitude is small. The 
mean value of £ from the data in iRedfield and Linsky I (120041 ) is 2.24 km/sec. The amplitude 
of turbulence has been discussed in Section 4.1.2 above, and for fluctuations of all kinds 
(parallel and perpendicular to a large scale field) can be defined as m = Sv/Va- We can use 
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the data in Table 1 to calculate the Alfven speed in the Local Clouds. The value obtained 
depends on whether the density used is the total density (ion plus neutral) or only the ionized 
component. The appropriate value depends on the scale of the fluctuations involved. For 
large-scale fluctuations (larger than the ion-neutral collisional mean free path) ion-neutral 
collisions cause the neutral atoms to be carried along in the Alfven wave, and the total density 
should be used. For small scale fluctuations (scales much smaller than the ion-neutral mean 
free path) the Alfvenic waves or fluctuations occur only in the ionized fluid, and the ionized 
density should be used. Since £ is an rms value formed from all fluctuations, the dominant 
contribution to £ should come from large scales of order the outer scale, and the relevant 
density for calculating the Alfven speed should be the total density. We will use the total 
density in the calculations below. 

Given the above considerations and the data in Table 1, the Alfven speed in the Local 
Clouds is 12.8 - 17.0 km/sec, with the range reflecting the possible values for the magnetic 
field strength. In calculating the dimensionless amplitude of the turbulence, Sv/Va, we 
assume Sv = \/3£, since £ corresponds to only one compon ent of the turbulent v elocit y 



fluctuations. For £ we use the mean of the entire sample of iRedfield and Linsky I (120041 ) . 
£ = 2.24 km/sec, giving 5v = 3.88 km/sec. We then have for the dimensionless amplitude 
of the turbulence 0.23 < Sv/Va — 0.30. Once again, for Alfvenic turbulence, Sb/B would 
have the same range. 

These calculations indicate that the turbulence in the Local Clouds is of substantial 
amplitude, but probably not so large as to cause significant excursions of the magnetic field 
direction, i.e. variations in field direction < 17°. These calculations argue that turbulence in 
these clouds is not sufficiently large to cause randomization of the field direction and thereby 
produce loss of intrinsic velocity anisotropy, if it were present. 



5.2. Limits on Ion Cyclotron Resonance Heating 

Perhaps the most striking difference between coronal and solar wind turbulence on one 
hand, and that in the Local Clouds on the other, is the absence in the Local Clouds of 
collisionless plasma processes responsible for enhanced heating of ions with larger Larmor 
radii (see Figure 6). The existence of a single temperature for many ions (and neutral atoms) 
is an observational result which can not be an "artifact" of a randomized interstellar magnetic 
field. In this section, we conjecture on the physical processes responsible fo r this apparent 



therm al equilibrium. A preliminary discussion of this topic has been given by lSpangler et al 



( 1201 If ) 
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The simplest way of explaining this result is to invoke collisionality. The solar corona 
is highly collisionless, and the solar wind at 1 au ranges from collisionless to only weakly 
collisional. Although the density of the Local Cloud plasmas is lower than either of the 
heliospheric plasmas we have discussed, the temperature is also lower, leading to higher 
collision frequencies. More importantly, the Local Cloud plasmas possess a significant neutral 
component (see Table 1), so ion-neutral collisions occur. The role of collisions is to take 
energy input by collisionless processes into one ion species, and into a limited number of 
degrees of freedom (i.e. perpendicular, but not parallel motions), and redistribute it to 
many species and all degrees of freedom. Observational support for the idea that co llisions 
will fulfill this role is given by solar wind observations reported by lKasper et al I ( 120091 ) . More 
collisional parts of the solar wind, such as the heliospheric current sheet, lack properties such 
as mass-proportional temperature and temperature anisotropy. We now consider whether 
the ion-neutral collision frequency in the Local Clouds is sufficiently large to make these 
plasmas collisional. 



Spangler et al I ( 120111 ) use two definitions of collisionality employed by lUzdensky I (120071 ) 



in the context of magnetic field reconnection in the solar corona. 



1. According to the first criterion, a collisionless plasma is one for which the ion cyclotron 
frequency is higher than the ion collision frequency. This permits, for example, insta- 
bilities with growth rates of the order of a fraction of the cyclotron frequency to develop 
without modification by collisions. 

2. In the second, much more restrictive criterion, a collisionless plasma is one for which 
the collisional mean free path is larger than the dimensions or characteristic scale of 
the plasma or medium. The converse situation of a collisional mean free path much 
smaller than the size of the system would then constitute a collisional plasma, even if 
the medium satisfied the first criterion for a collisionless plasma given above. 



To evaluate the co llisionality of th e Loca l Cloud plasmas, we need to identify the relevant 
collisional processes. ISpangler et al I (1201 if ) discuss charge exchange and induced dipole 
scattering. These two processes have similar collision cross sections, with charge exchange 
being larger by a factor of a few for Local Cloud conditions. The microphysics of what 
occurs is quite differe nt in the two cases. For conditions appropriate to the Local Clouds, 



Spangler et al I ( 120111 ) calculate an H + -H ion-neutral collision frequency due to induced dipole 
scattering of 3 x 10 -10 Hz. This frequency i s about 8 orders of m agnitude smaller than the 
proton ion gyrofrequency (see Table 2 of ISpangler et al 1 120111 ). meaning that the Local 
Clouds easily satisfy the first of the above criteria for being collisionless. However, the 
collisional mean free path corresponding to this collision frequency is 5 x 10 15 cm = 1.5 x 10 -3 
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parsec = 330 AU. For the case of H + -H charge exchange, the collision frequency would 
be a few times higher, and the collisional mean free path a few times smaller than the 
numbers given immediately above. In any case, we conclude that the Local Clouds are 
highly collisional by Uzdensky's second criterion. 

It seems likely that collisions between ions taking part in the Alfvenic waves and tur- 
bulence in the local clouds, and neutral atoms which are not, is responsible for the removal 
of temperature anisotropy and ion-specific heating. It would be worthwhile to conduct a 
theoretical study of ion and neutral atom dynamics in the presence of Alfvenic turbulence, 
including collisions due to both charge exchange and induced dipole scattering. Such cal- 
culations could determine if the collisions do indeed remove temperature anisotropy and 
ion-specific temperatures, by taking the energy which is preferentially input to one or a few 
ions, and distributing it to all ion and neutral atom species. 

Introducing Uzdensky's first criterion of collisionality has provided interesting insight 
in the present context. The fact that the collision frequency is approximately 8 orders of 
magnitude smaller than the gyrofrequency means that plasma instabilities and quasilinear 
modification of ion distribution functions would have ample time to develop, unencumbered 
by collisions, if conditions for instabilities were present. The absence of "collisionless" obser- 
vational signatures suggests that energetic processes on ion temporal and spatial scales are 
not occurring in the Local Clouds. 

6. Summary and Conclusions 

The conclusions of this paper are as follows. 

1. High-resolution optical and UV spectroscopy of absorption lines for several ions and 
atoms in the Local Clouds of the VLISM provide a remarkable number of diagnostics 
of the turbulence in these clouds. The data are sufficient to make some tests of the 
similarities to, or differences from, MHD turbulence in the solar corona and solar wind. 

2. The turbulent line broadening parameter £ does not show a systematic dependence on 
direction in the sky, as would be expected if the velocity fluctuations are predominantly 
transverse to the local Galactic magnetic field. This result indicates either that the 
turbulent fluctuations are not anisotropic with respect to the local magnetic field, or 
that the field is partially randomized on a size scale of several parsecs. 

3. The quantitative limits we can place on anisotropy in the case of an ordered field are, 
however, not very restrictive. Our results indicate an upper limit to the anisotropy 
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parameter e < 0.70. If such a large value of the anisotropy is present, it must be 
masked by another, unidentified process that causes variations in £ from one line of 
sight to another, which in many cases are within the same cloud. 

Although a conclusion of this paper is that there is no significant anisotropy in £, 
and therefore no result for the direction of the local interstellar magnetic field Bo, a 
set of candidate directions for B$ c hosen for further analy sis are in rough agreement 



with the field direction proposed by lLallement et al I ( 120051 ) . This could justify future 



investigation with a larger set of absorption line measurements. 

The temperature T obtained from the spectral line widths (see Equation (1)) also 
shows no dependence on direction on the sky. The most straightforward interpretation 
of this result is that there is not a temperature anisotropy T± > T\\ as is the case in 
the solar corona, and to a lesser extent, the solar wind. 

We determined an upper limit to the amount of ion-mass-dependent heating by study- 
ing the statistical acceptability of Equation (10), a temperature model with mass- 
dependent temperature. This model was fit to the data for 9 lines of sight with 7 or 8 
transitions from ions or atoms of different mass. The maximum value of d allowed for 
these 9 lines of sight was d = 0.74, with the other lines of sight having smaller values 
of d. The case of mass-proportional temperature would have d = 1. 

The explanation for results # 3-5 is most likely that the turbulence (and ion and neutral 
atom distribution functions) have persisted for many ion-neutral collision timescales. 
These ion-neutral collisions have apparently eliminated temperature anisotropics and 
temperature differences between species. This conclusion itself is of interest in that 
it indicates the absence of energetic kinetic plasma physics processes on the scale of 
the ion cyclotron radius or ion inertial length in the Local Clouds of the Very Local 
Interstellar Medium. 



This work was supported at the University of Iowa by grants AST09-07911 and ATM09- 
56901 from the National Science Foundation. 



A. A model for turbulent fluctuations 

In this appendix, we derive the relationship between the probability distribution function 
of velocity fluctuations in a frame of reference defined by the interstellar magnetic field, and 
the distribution function in an observer-centered frame defined by the line-of-sight to a star. 
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The anisotropic turbulent velocity fluctuations are most naturally defined in a coordinate 
system in which one coordinate, the z coordinate, is in the direction of the large scale 
magnetic field B . Let b be the unit vector which points in the direction of B . The unit 
vector b can be defined by the galactic longitude and latitude (A, (3) towards which it points. 
In the Bq coordinate system, the z axis is in the direction of b, x is in the plane defined by 
b and the direction to the north galactic pole, and y completes a right-handed coordinate 
system. 

The turbulence model we want to test is one in which the fluctuations in the plane 
perpendicular to b are larger than those in the direction of b. A probability distribution 
function which is simple in mathematical form and describes this is 

1 \ , , vl + vl 



™ = Kwwwi) exp( "if )exp( -^vf ) (A1) 

with V± > V\\ by assumption. The distribution function (Al) satisfies the normalization 
requirement J d 3 vf(v) = 1. 

The observed spectral line width is proportional to the rms fluctuation in the component 
of the velocity along the line of sight, which is in the direction of the unit vector I. The unit 
vector / points in the direction of galactic longitude and latitude (l,b). In what follows, we 
assume that the observed turbulent width £ can be expressed as 

£ 2 =< v\ >=< v 2 L (v x , v y , v z ) >= J d 3 vv 2 L (v x , v y , v z )f(v) (A2) 

where vl is the line-of-sight component of the flow velocity. By expressing £ 2 as an expec- 
tation value, we assume that the line of sight integration through the cloud samples a large 
number of independent eddies in the cloud and thus satisfies the ergodic theorem. 

To evaluate the expression (A2), we need to express vl(v x , v y , v z ). We define a coordinate 
system such that the unit vector for one coordinate is in the direction of I. The unit vector 
corresponding to another coordinate (e^) is in the plane defined by and the direction to 
the north galactic pole. Finally, the unit vector for the third coordinate in the line-of-sight 
coordinate system, ip, is defined by ct x ep = e^. 

The transformation between the two coordinate systems is given by 

(A3) 







(v, 


Vp 


|.T.| 




vl J 







where T is the matrix which generates the rotation from the magnetic-field oriented coordi- 
nate system to the line-of-sight coordinate system. The matrix T is given by an Euler angle 
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transformation, so it is the product of three matrices, each describing the rotation through 
one of the Euler angles, 

T = T^®T 2 ®T 1 (A4) 

The Ti operator is the first one. It rotates the (x, y, z) coordinate system about the y axis 
so that the z axis is in the galactic plane. It is given by 

/ cos/3 sin/3 \ 
Ti=l 10 (A5) 
\ - sin/3 cos/3 / 

The next rotation operator T 2 rotates the coordinate system in the galactic plane such that 
the z" axis is pointing along the galactic longitude / of the line of sight. It effects a rotation 
about the x axis through an angle Al = A — I 

1 
/ | cos Al -sin Al ) (A6) 
sin Al cos Al 

The third rotation operator rotates the coordinate system "upwards" about the y" axis 
through an angle b. At this point, the z" direction coincides with I. The operator T 3 is given 
by 

(cos b — sin b \ 
10 (A7) 
sin b cos b J 

Multiplying the three matrices together as in Equation (A4) to produce T, we have 

(cos b cos j3 + cos Al sin b sin (3 — sin Al sin b cos b sin [3 — cos Al sin b cos (3 \ 
sin Al sin (3 cos Al -sin Al cos (3 (A8) 

sin b cos (3 — cos Al cos b sin (3 sin A/ cos 6 sin b sin /3 + cos A/ cos 6 cos (3 J 

We are interested in only the vl component in the (T,P,L) coordinate system, which is given 
by 

v L = (sin b cos (3 — cos Al cos b sin /3)v x + (sin A/ cos b)v y + (sin 6 sin (3 + cos A/ cos b cos /3)t> 2 

(A9) 

In keeping with our model of the clouds and their turbulence, Equation (A9) holds at every 
point along the line of sight through a cloud. The quantity which is measured as the turbu- 
lence parameter, £ 2 =< v\ >. We assume that < v x v y >=< v x v z >=< v y v z >= (i.e. the 
fluctuations have no dominant polarization). We then have 

£ 2 =< v\ >= (sin b cos (3 - cos Al cos b sin (3) 2 < v\ > +(sin Al cos b) 2 < v 2 > (A10) 

+ (sin b sin f3 + cos Al cos b cos (3) 2 < v 2 > 
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The expectation values of the squares of the velocity components are 

< v 2 >= J d*vvff(v) (All) 

% = x,y,z, and for our model < vl >=< v 2 >. Substituting Equation (Al) into (All) and 
evaluating gives 

< vl >=< vl >= V[ (A12) 
< v 2 z >= If (A13) 

v 2 

We now define a first anisotropy parameter 77 = and use it in Equation (A10) to give 

£ 2 < v 2 > 

— 2 = — -j^ — = [(sin b cos f3 — cos Al cos b sin j3) 2 + (sin Al cos b) 2 + ^(sin 6 sin (5 + cos A/ cos b cos /3) 2 ] 

(A14) 

This is an appealing equation. The measured square of the turbulent velocity width is on 
the left hand side of the equation, and functions of galactic coordinates (of the line of sight 
and the large scale B field), as well as the anisotropy parameter 77 are on the right. However, 
the expression can be simplified further by use of trigonometric identities, and definition of 
a second anisotropy parameter e, r\ = 1 — e, to be 

e <v\> 



V 2 ~ V 2 



= 1 - e(sin b sin (3 + cos Al cos b cos (5) (A15) 



This is one of the two fundamental equations in our analysis. Equation (A15) relates the 
observed turbulent line width £ (left hand side, normalized by the unknown but estimatable 
Vl) to the anisotropy parameter e, the direction of the line of sight (/, b) and the direction 
of the local interstellar magnetic field (A, (5) on the right hand side. Although this form of 
the equation is most useful for the analyses carried out in this paper, it is more instructive 
to rewrite Equation (A15) in terms of an angle A, defined by I • b — cos A. It is easy to show 
that Equation (A15) becomes 

f 2 < v 2 > 

|j = -^ = l-"os 2 A (A16) 

One of the goals of this paper is to determine if, for some set of values of (A, /3), Equation 
(A15) provides a good representation of the measured values of £ for some value of e > 0. 
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